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We consider the A'^-components 3-dimensional massive Gross-Neveu model compactified in one 
spatial direction, the system being constrained to a slab of thickness L. We derive a closed formula 
for the renormalized L-dependent four-point function at vanishing external momenta in the large- A*' 
limit (the effective coupling constant), using bag- model boundary conditions. For values of the fixed 
coupling constant in absence of boundaries A > Ac — 19.16, we obtain small-distance asymptotic 
freedom (for L ^ 0) and a singularity for a length L''"^ such that 2.07 < L'''' < 2.82 m~^, 
m being the fermionic mass. Taking for m an average of the masses of the quarks composing the 
proton, we obtain a "confining" length Lp"^' which is comparable with an estimated proton diameter. 

PACS number(s): ll.lO.Jj; ll.lO.Kk; 11.15.Pg 
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I. INTRODUCTION 



The currently accepted theory of strong interactions. 
Quantum Chromodynamics (QCD), has an intricate 
structure which makes it difficult to be directly handled. 
For this reason, simplified models have been largely em- 
ployed over the last decades in trials to obtain clues about 
the behavior of interacting quarks. Among them one of 
the most celebrated is the Gross-Neveu model Q , which 
considers the direct four-fermions interaction. In this pa- 
per, we consider the iV-component tridimensional mas- 
sive Gross-Neveu model, in the large- TV limit, compacti- 
fied along one of the spatial directions. From a physical 
point of view, the model is intended to describe fermions 
bounded between two parallel planes, a distance L apart 
from one another. We will call from now on quarks the 
fermions associated to the fermionic field, being under- 
stood that we are dealing with the quanta resulting from 
the quantization of the Gross-Neveu model. We define 
from the four-point function, an effective coupling con- 
stant g{3, L ; A) in the large- limit, which presents dif- 
ferent behaviours if the fixed coupling constant A in un- 
bounded space is below or above some "critical" value Ac. 
For A < Ac (weak coupling) the model presents small- 
distance asymptotic freedom, g{i, L ; A) ^ 0, as L ^ . 
For strong coupling (A > Ac), starting from small values 
of L, we analytically demonstrate the simultaneous exis- 
tence of short distance asymptotic freedom and of a sin- 
gularity for L equal to some L = L^'^\ g{3, L ; A) — > cx), 

(c) . 

as L L_ . This means that in the strong coupling 
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regime, if we start with a system of a quark- antiquark 
pair bounded between two planes a distance L < L^'^\ 
from one another it would not be possible to separate 
them a distance larger than L'^'^\ This spatial confine- 
ment of the quark- antiquark pair could be interpreted as 
the existence of bound states, characteristic of the model 
in the strong coupling regime. We will refer later to this 
property of the model as confinement, understood as the 
spatial confinement described above, not of colour con- 
finement as it should happpens for QCD. An interesting 
point is that the confining length obtained is comparable 
with an estimated classical proton diameter. 

Although perturbatively non renormalizable for di- 
mensions D > 2, the massive Gross-Neveu model in 
tridimensional space has been shown to exist and has 
been constructed in Ref. . The idea that perturbative 
renormalizability should not be an absolute criterion to a 
model to be physically consistent comes from works done 
more than one decade ago. The results in renormalization 
theory reported in Refs. 0,0] give a theoretical founda- 
tion to support the idea that non perturbatively renor- 
malizable models can exist and have a physical mean- 
ing 5j. In fact, the Gross-Neveu model for dimensions 
greater than D = 2 has been investigated since the mid- 
dle of the seventies. Since then we know that the diffi- 
culties related to the perturbative non-renormalizability 
of the model can be surmounted, at least for D = 3, by 
performing partial summations in the bare perturbation 
series^while maintaining its initial structure (see Refs. 
QiIEQ)' This point has been also studied in Ref. 0. 
The essence of the method lies in the fact that the model 
in three dimensions may acquire a naive power counting 
behavior of a just renormalizable theory, if we perform 
a partial summation, to all orders in the coupling con- 
stant, of the chains of one- loop bubbles before trying to 
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renormalize the perturbation series directly. 

In present work, we perform such a summation for the 
large- model and we renormalize it employing a modi- 
fied minimal subtraction scheme, where the coupling con- 
stant counter-term is a pole at the physical value of the 
argument of generalized Epstein zeta-functions 0. This 
technique, which employs analytic and zeta-function reg- 
ularizations, has been recently used to study the bound- 
ary dependence of the coupling constant and the sponta- 
neous symmetry breaking in the compactified 0* model 

caini. 



constraining the i-dependence of the field to a segment 
of length L, allow us to proceed with respect to the x 
-coordinate in a manner analogous as it is done in the 
imaginary-time formalism in field theory. This extension 
of the Matsubara procedure to describe a model compact- 
ified in a space-coordinate has been already used in the 
case of bosonic fields ^3 ^3 ■ Accordingly, for fermionic 
fields, the Feynman rules should be modified following 
the prescription 
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II. THE COMPACTIFIED MODEL 

We start from the Wick-ordered massive Gross-Neveu 
Lagrangian in Euclidian metrics (in relativistic units, h — 
c=l), 

£ =: ^(x)(*7''a^ + m)V'(a^) : +\ ■ ' :, (1) 

where m is the mass, u the coupling constant and a; is a 
point of . The fermionic field of spin i has N compo- 
nents, 'ip'^(x) a = l,2,...,iV, and summation over fiavor 
and spin indices are understood. For a finite number of 
flavors A^, for 2 < D < 4, a throughout study has been 
done in Ref. Here we consider the limit N ^ oo, 

which permits considerable_simplifications. 

In the language of Ref. 2] a theorem has been proved 
stating that a family of parameters can be exhibited such 
that, for iV sufficiently large, all renormalized Schwinger 
functions of the type 52p(m, A, iV), p=l,2,... exist. Thus 
we examine in the following the large- iV zero external- 
momenta four-point function, which is the basic object 
for our definition of the renormalized coupling constant. 
The four-point function at leading order in is given by 
the sum of all chains of single one-loop diagrams. This 
sum gives for the four-point function at zero external 
momenta the formal expression Q 



i + Nuj:{d)' 
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where S(£)) is the Feynman integral corresponding to the 
single one-loop four-point subdiagram 
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Notice that has dimension of raass^^"^ and the 

coupling constant u has dimension of mass^^^ . 

We consider the system bounded between two parallel 
planes, normal to the x -axis, a distance L apart from one 
another. We use Cartesian coordinates r — {x,z), where 
z is a {D— l)-dimensional vector, with corresponding mo- 
menta k — {kx, q), q being a {D — l)-dimensional vector 
in momenta space. In addition, we assume that the field 
satisfies bag-model boundary conditions 0, on the 
planes x = and x = L. These boundary conditions. 



Then, the L-dependent four-point function at zero exter- 
nal momenta has the formal expression. 



rg)(o,L) = 



1 + iViiS(D,L)' 



(5) 



where S(-D, L) is the L-dependent Feynman sum-integral 
corresponding to the single one-loop subdiagram. 
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Using the dimensionless quantities qj = kj 127:7X1 (j — 
1,2, ...,D-1) anda = {mL)~^, the single bubble T,{D,L) 
can be written as 



S(L>,a) 



where 
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U3{s,a) = 
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These expressions were obtained by performing the in- 
tegrals over Qj in Eq. ^ using well-known dimensional 
regularization formulas {fj] . 

Sums of the type appearing in Eq. © can be casted 
in the form 
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where 
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is one of the generalized Epstein zeta-functions de- 
fined for Jiev > 1/2, which possesses the following ana- 
lytic extension 
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valid for the whole complex i^-plane. Therefore, the func- 
tions Ui{s,a) {I = 1,2,3) can be also extended to the 
whole complex s-plane, and we obtain 
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where 
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Replacing these expressions in Eq.Q, we see that 
I](I?, L, s) can be written in the form 



£(£>, L, s) = H{D.s) + G{D, L, s) 



(19) 



where H{D, s) is a term independent of a, coming from 
the first parcels between the brackets in Eas. p4|l . p 51) and 
(|16|l . while G{D,L,s) is the term arising from the sum 
over the parcels containing the W^-functions in these 
equations. We are using a modified minimal subtrac- 
tion scheme, where the coupling constant counter-terms 
are poles appearing at the physical value s = 2. Then 



the L -dependent correction to the coupling constant is 
proportional to the regular part of the analytical exten- 
sion of the Epstein zeta -function in the neighborhood 
of the pole at s = 2. We see from Eqs. ((T^ and 
that, for s = 2 and even dimensions D > 2, H{D) is 
divergent due to the pole of the F-function. Accordingly 
this term should be subtracted to give the renormalized 
single bubble function Ei^(Z), L). For sake of uniformity, 
the term H{D) is also subtracted in the case of other di- 
mensions D, where no poles of F-functions are present; in 
these cases, we perform a finite renormalization. Then, 
we simply get, for any dimension D, 
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In dimension D — 3, using the explicit form for the 
Bessel functions Kj-i/2{z) = y/ne~^ /\/2z and perform- 
ing the resulting sums, we obtain (remembering that 
a = {mL)~'^) 
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Substituting this expression into Eq. (jSJ, we obtain the 
large- effective renormalized coupling constant, taking 
Nu — A fixed as usual, 



g{3,L)^NTi'^{0,L) 
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III. BOUNDARY BEHAVIOUR 

Fig. 1 shows the renormalized single bubble (in units 
of m), S — I]fl,(3, L)/m, as a function of L (in units 
of m~^). We see from this plot and from Eq. H21() 
that lim 00^ r{S, L) = and that I]_r(3, L) diverges 
as L ^ 0. This means that (see Eq. (|22|) ') the fixed 
coupling constant A corresponds to the large-A^ renor- 
malized effective coupling in absence of boundaries and 
that, independently of the value of A, the system presents 
ultra-violet asymptotic freedom. Moreover, E/j(3, L) has 
a zero at i = -^mL — 2.07, for which g{3, L) = A. 

The interesting point is to investigate what happens 
for finite values of L, when we start from L near zero 
(in the region of asymptotic freedom). In Fig. 2, we have 
plotted the relative coupling constant y = g{3,L)/X as 
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FIG. 1: Renormalized single bubble (in units of m), S 
T,r{3, L)/m, as a function of the length L (in units of 



a function of L (in units of to~^), taking four different 
values of A (also in units of m~^). For the smallest value 
considered (A — 0.2), we only see asymptotic freedom 
with g{2>,L) varying from (for L = 0) up to A (for L — 
oo). Increasing the value of A (e.g. A = 10 and A = 15 in 
Fig. 2), all the curves present a peak at L = Lmax — 2.82, 
the point at which X/g{3, L) reaches its minimum value. 
This minimum value of X/g{3,L) vanishes for A = Ac — 
19.16 (corresponding to the full line in Fig. 2) and so 

the effective coupling constant g(3, Lmax) diverges. This 
means that, for this value of A, the system is absolutely 

(c) 

confined in the slab of thickness Lmax- Taking greater 
values of A > Ac, the minimum of X/g(3,L) becomes 

negative, g{3, L) goes to oo at a value L^'^^ < Lmax and 
the system would be confined in a thinner slab. However, 
as A goes to oo, L'^'^^ tends to a lower bound equal to 

(c) 

L^-„. In other words, for Ac < A < oo, the effective 
coupling constant is singular for a critical value of L, 
namely L*^'^)(A); the system is then confined in a slab of 



thickness L('^'(A) lying in the interval ^i|^)„,i 



max 



Searching for a physical interpretation of our results, 
we shall take the fixed coupling constant A = Ac imply- 
ing that the physical confining length is L'^'^^ = Lmax, 
the maximum conceived confining length. In order to 
estimate the value of L^'^^ we need to fix a value of 
TO. We choose to = 3.25 x IO^^cto"^ 6.5Meu), 
corresponding to an average between typical u- and d- 
quark masses (recall that the quark content of the pro- 
ton is {uud}) |l0i. This choice gives a confining length 
L^p'^ ~ 8.67 X 10-12 cm. This value should be com- 
pared with an estimated classical proton diameter, 2Rp ^ 
6.92 X \Q~^'^cm notice that this value is about lO'^ 
times larger than the experimentally measured proton 
charge radius jl3- 



FIG. 2: Plots of the large-A'^ relative effective coupling con- 
stant, y = (ji(3, L)/\, as a function of L (in units of m~^) for 
four different values of the fixed coupling constant A (in units 
of m~^): A = 0.2 (dashed line), A = 10 (dotted line), A = 15 
(dotted-dashed line) and A = 19.16 (full line). The dotted 



vertical lines, passing by 1/^^^ 
plotted as a visual guide. 



2.07 and L 



2.82, are 



IV. CONCLUDING REMARKS 



Under the assumption that the large-iV Gross-Neveu 
model can be thought about as a model sharing some 
of the properties of physically meaningful theories (in 
the sense of experimentally testable), we can say that we 
have exhibited in a direct way the properties of ultravio- 
let asymptotic freedom and confinement, characteristics 
of strong interactions. This could be contrasted with the 
fact that in QCD, although asymptotic freedom is a well 
established property on theoretical grounds, confinement 
has been confirmed only by computer simulations in lat- 
tice field theory. Confinement has never been, to our 
knowledge, analytically proved. 

It is worth mentioning that, in Ref. p^. we have an- 
alyzed the dependence of the coupling constant on L for 
the large-A^ (j)'^ model, the bosonic counterpart of the 
Gross-Neveu model studied here. In that case of bosons, 
no confinement is found at all, but only asymptotic free- 
dom is encountered. Comparison of the results derived in 
Ref. [1.0] with those of the present work suggests that the 
effect of confinement in hadronic matter may be mainly 
due to the fermionic nature of the quarks. 

Finally, we would like to point out that the behavior 
described in this work should not be heavily dependent 
on the confining geometry. One expects to have ultravio- 
let asymptotic freedom in any case and partial or absolute 
confinement, depending on the value of A. A thorough 
study on the subject with different geometries, and in- 
cluding temperature, is in progress and will be presented 
elsewhere. 
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